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FOREWORD 


Work  on  the  characterization  and  mathematical  modeling  of  classes  of 
PCM  Bit  Synchronizers  was  conducted  under  AIRTASK  A5355352-054E- 
3F09905003  to  provide  technical  support  to  the  Telemetry  Group  of  the  RCC 
(Range  Commanders  Council).  The  work  is  part  of  Work  Unit  Assignment 
A5352 10000002,  Range  Telemetry  Support.  The  purpose  of  the  task  is  to 
conduct  a  program  of  testing  and  investigation  that  will  provide  information 
which  can  be  used  by  the  Telemetry  Group  to  keep  the  RCC  publications  of 
"Telemetry  Standards"  and  "Test  Methods  for  Telemetry  Systems  and  Sub¬ 
systems"  updated  and  abreast  of  technological  advances. 

The  investigations  were  conducted  and  the  report  prepared  by  Dr.  William 
C.  Lindsey,  University  of  Southern  California,  under  Navy  contract  N61756- 
72-M-6443. 
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SUMMARY 


This  report  presents  the  results  obtained  during  the  study  on  mathe¬ 
matical  modeling  of  PCM  bit  synchronizers  {also  known  as  PCM  signal 
conditioners).  It  is  written  in  two  parts.  The  first  part  presents  results 
obtained  on  the  study  of  the  effects  which  bit  synchronization  has  on  the  bit 
error  rate  (BER).  Several  PCM  baseband  signaling  formats  are  included 
in  the  study,  viz.,  RZ,  NRZ,  Manchester,  and  Miller  Coding.  Theoretical 
performance  limitations  of  bit  synchronizers  and  the  associate  bit  detector 
are  developed  by  relating  the  BER  to  the  signal-to-noise  ratio  per  bit.  It 
is  assumed  that  the  noise  added  to  the  input  signal  is  white  and  Gaussian, 
the  PCM  signal  waveforms  are  either  rectangular  or  filtered  and  that  bit 
synchronization  is  perfect.  Next  the  constraint  of  perfect  bit  sync  is 
relaxed  to  include  the  more  practical  case  where  bit  sync  jitter  and  bit 
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slippage  affect  the  detection  of  the  data  bite.  To  account  for  the  effects 
of  bit  sync  jitter  on  the  BER,  the  probability  density  function  of  bit  sync 
jitter,  valid  for  a  broad  class  of  bit  synchronisation  systems  (BSS),  is 
postulated  and  this  is  interconnected  with  the  theory  of  bit  detection.  The 
deleterious  effects  are  described  in  terms  of  the  parameter  o  =  To^; 
where  a ^  represents  the  normalized  bit  sync  jitter  and  T  is  bit  duration 
in  seconds.  Various  design  curves  and  performance  comparisons  are 
made  for  several  PCM  baseband  signaling  formats. 

The  second  part  of  this  report  is  devoted  to  the  study  of  bit  synchroni¬ 
zation  performance  characterization.  A  general  discussion  is  given  which 
highlights  the  many  performance  measures  which  serve  to  characterize  the 
performance  of  BSS.  Some  general  mathematical  relations  are  given  for 
evaluation  of  the  acquisition  time,  acquisition  range,  and  bit  slippage  rate 
for  PCM  bit  synchronizers  which  are  mechanized  sucn  that  the  phase 
detector  characteristic  is  periodic. 

This  theory  is  followed  by  a  discussion  of  bit  synchronizer  performance 
comparisons  and  this  leads  to  the  suggestion  of  the  mathematical  models 
which  can  be  implemented  using  the  method  of  computer  software.  The 
four  parameters  of  interest  in  the  computer  software  model  include  bit 
sync  acquisition  time,  bit  slippage  rate,  bit  sync  jitter  and  bit  sync 
acquisition  range.  Two  bit  sync  approaches  which  will  lead  to  a  specific 
software  package  are  recommended:  These  include  the  absolute  value, 
least  squares  and  decisions  directed  approach.  The  approach  leads  to 
suggestions  of  techniques  which  can  be  used  in  the  procurement  of  BSS 
as  well  as  the  development  of  standard  procedures  for  testing  PCM  BSS« 
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I.  THEORETICAL  PERFORMANCE  OF  IDEAL  PCH  SIGNAL  CONDITIONERS* 


In  this  part  we  consider  the  problem  of  detecting  two  equally  probable 
signal  waveforms  s  (t)  or  s  (t),  each  of  duration  T  seconds,  in  a  back- 
ground  of  white  Gaussian  noise  of  single-sided  spectral  density  of 
watts /Hz.  The  various  PCM  signaling  formats  to  be  considered  are 
illustrated  in  Fig.  1-1. 

Now  suppose  that  in  each  of  a  succession  of  time  intervals  of  duration 
T  we  know  that  the  received  signal  consists  of  either  signal  s^(t)  or  s^t) 
and  additive  noise,  say  n(t).  We  seek  a  data  detector  (wherein  bit  sync 
is  assumed  known)  that  will  decide  on  the  basis  of  the  received  waveform 
in  each  interval  which  of  the  two  possible  signals  was  transmitted.  If 
we  choose  the  figure  of  merit  for  such  a  device  to  be  the  error  probability 
then  the  optimum  (perfectly  synchronized)  data  detector  is  one  that  will 
announce  s^(t)  if 

T  rp 

Tf  f  [x(t)  -  s  (t)]2dt<^  f  [x(t)  -  s  (t)]2dt  (I- 1 ) 

T  J0  1  Tv0  & 

Here  x(t)  represents  the  received  signal  plus  noise.  Assuming  s^(t)  and 
s^ft)  have  equal  energy  i.  e.  , 

T  T 

E  =  f  s  (t)dt  =  I  s  (t)dt  =  ST  (1-2) 

Jo  1  Jo  2 


*In  this  report  the  term  "signal  conditioner"  is  used  to  imply  the  bit 
detector  and  bit  synchronizer  as  a  combined  operational  system. 


3 


1 


1 


1 


1 


I - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 

I  1  I  0  •  1  •  1  •  0  •  0  1  0  1  1  1  1  1  0  1  1  • 

I  I  I  I  I  l  l  l  l  l  l  l 


I  I 


I  I 


|  NRZ-LEVEL  (OR  NRZ  CHANGE). 

"ONE"  IS  REPRESENTED  BY  ONE  LEVEL. 


I- 


ZERO"  IS  REPRESENTED  BY  THE  OTHER 


LEVEL. 


I  I 


I  I  I 


I  I 

I  I  I 


I  I 


I  I 

I  ! 

i  i 

I  I  I 


I  I 


I  I  i  I 


1  1  1 

i  i  i 

1  r 

u  i  U 

i  i  l 

i  i  i 

U  i  i  i 

lil 

j  NRZ-MARK  (DIFFERENTIAL  ENCODING). 

"ONE"  IS  REPRESENTED  BY  A  CHANGE 
j  IN  LEVEL. 

j  "ZERO"  IS  REPRESENTED  BY  NO  CHANGE 
|  IN  LEVEL 

I 

I 

-NRZ -SPACE  (DIFFERENTIAL  ENCODING) 

I  "ONE  ’  IS  REPRESENTED  BY  NO  CHANGE 
|  IN  LEVEL. 

|  "ZERO"  IS  REPRESENTED  BY  A  CHANGE 
|  IN  LEVEL. 

I 


I 


U1 


I  I  I  I 
I  I  I  I 
1  I  I  I 


I  I 


I  I  I 


I  I 


RZ. 

"ONE"  IS  REPRESENTED  BY  A  HALF-BIT 
I  WIDE  PULSE. 

|  "ZERO"  IS  REPRESENTED  BY  NO  PULSE 
|  CONDITION. 


I  I  I  I  I 


unmjnjimjimiji 


I  I  I 

i  i  i 

i  *  i  i 


i 


■  BI  PHASE-LEVEL  VSPLIT  PHASE  OR 
(MANCHESTER  CODE). 

"ONE"  IS  REPRESENTED  BY  A  HALF-BIT 
I  WIDE  PULSE  OF  ONE  POLARITY. 

I  "ZERO"  IS  REPRESENTED  BY  A  HALF-BIT 
WIDE  PULSE  OF  THE  OPPOSITE  POLARITY. 


iiuinrmjuLnj 


i  i 


i  i 


i  i 


i  i 

i  i 

i  i 


uirumrurjin 


i  i  t 

i  i  i 

i  i  i 

i  i  i 

!  I  T 

I  I  I 

U  I  ! 


i  i 


I  I 


i  i  i 
i  i  i 
i  i  i 
i  i  i 


BI-PHASE  MARK. 

A  TRANSITION  OCCURS  AT  THE  BEGINNING 
OF  EVERY  BIT  PERIOD. 

"ONE"  IS  REPRESENTED  BY  A  SECOND 

|  |  TRANSITION  ONE-HALF  BIT  PERIOD  LATER. 

j  "ZERO"  IS  REPRESENTED  BY  NO  SECOND 
(TRANSITION. 

I  I 

i  i 

Bl -PHASE  SPACE. 

|  A  TRANSITION  OCCURS  AT  THE  BEGINNING 
OF  EVERY  BIT  PERIOD. 

"ONE"  IS  REPRESENTED  BY  NO  SECOND 
!  TRANSITION. 

'  "ZERO"  IS  REPRESENTED  BY  A  SECOND 
j  TRANSITION  ONE-HALF  BIT  PERIOD  LATER. 

|  DEI  AY  MODULATION  (MILLER  CODE). 

"ONE"  IS  REPRESENTED  BY  A  SIGNAL  TRAN- 
I  SITION  AT  THE  MIDPOINT  OF  THE  BIT  PERIOD. 
"ZERO"  IS  REPRESENTED  BY  NO  TRANSITION 
UNLESS  IT  IS  FOLLOWED  BY  ANOTHER  ZERO. 


1  I 

i 

L 


IN  THE  LATTER  INSTANCE.  A  TRANSITION  IS 
[  PLACED  AT  THE  END  OF  THE  BIT  PERIOD  'jF 
I  THE  FIRSI  ZERO. 


Figure  l-l.  I ’CM  Signaling  Formats. 
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This  optimum  decision  rule  can  be  realized  (see  Fig.  1-2)  by  constructing 
two  cross-correlators  and  providing  perfect  bit  synchronization  to  be 
used  to  sample  the  cross-correlator  outputs.  If  the  first  cross-correlator 
output  exceeds  the  second,  we  decide  that  s^(t)  was  present;  otherwise, 
decide  s^t).  Alternatively,  the  optimum  decision  circuit  could  be  realized 
by  constructing  two  matched  filters  having  impulse  responses  Sj(T-t) 
and  s^(T-t);  however,  we  do  not  pursue  this  realization  since  one  generally 
implements  Fig.  1-2  in  practice. 


Figure  1-2.  Optimum  Decision  Circuit  Assuming  Perfect  Hit  Syne. 


I-  1.  Theoretical  Performance  of  Perfectly  Synchronized  PCM  Signal 
Conditioners.  The  measure  of  system  performance  is  bit  error  rate 
(BER).  This  measure  of  performance  represents  the  average  number  of 
times  an  incorrect  bit  state  occurs  in  a  given  time  period.  Frequently, 
the  BER  is  expressed  as  the  bit  error  probability  (BEP)  which  is  the 
ratio  of  the  number  of  bits  in  error  to  a  given  number  of  bits  transmitted. 
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In  what  follows  we  shall  denote  the  BEP  by  the  symbol  Pg.  For  the 
case  of  perfect  bit  synchronization  the  BEP  turns  out  to  be  a  function  of 
the  ratio  of  the  signal  energy  per  bit,  E  =  ST,  to  single- sided  noise  spectral 
density  =  kT^  watts/Hzl  This  ratio,  say  R  =  E/Nq  =  ST/N^,  will  be 
referred  to  as  the  signal  energy-to-noisc  ratio  per  bit  or  more  simple 
the  signal-to-noise  ratio  (SNR).  We  now  present  the  theoretical  perform¬ 
ance  limitations  for  various  PCM  signal  conditioners  relating  BEP  to  the 
signal-to-noise  ratio  R  =  E/N^  =  ST/N^. 


(a)  The  RZ  Signaling  Format.  In  what  follows  we  shall  summarize 
the  performance  of  two  types  of  RZ  signaling.  In  both  cases  the  signals 
Sj  and  s^  are  selected  such  that  they  are  orthogonal,  i.  e.  , 


J  s1(t)s2(t)dt  =  0 

Type  I-RZ  signals  are  represented  by 

s  (t)  =  0  s  t  £  T 

1 


s  (t)  =  0 


(1-4) 


(1-5) 


The  BEP  for  Type  I  -  RZ  signaling  is  easily  shown  to  be  given  by 


PE  =  j  erfc(v^74) 


(1-6) 


*The  thermal  noise  power  generated  per  unit  bandwidth  is  v2  =  4kT°R.  When 
this_source  is  matched  to  an  ideal  load  of  R  ohms  the  noise  power  delivered 

is  v2/(R/2R)2. 
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where 


A  2  r 

erfc(x)  -  ^-r  | 


.t2 

e  dt 


(1-7) 


This  assumes  that  a  threshold  is  established  in  the  data  detector  at  one- 
half  the  signal  energy.  If  >  ST/2  then  s^  is  announced;  otherwise,  a 
decision  is  made  in  favcr  of  s?.  This  threshold  is  optimal  for  correlation 
detection. 


Type  II  -  RZ  signals  are  defined  by 

0st^  T/2 


s  jft)  = 


-VS*  T/2  M^T 


s  <t)  =  0  *  t  *  T 


(1-8) 


while  the  corresponding  BEP  is  given  by 


P_  =  1  erfc(^72) 


(1-9) 


From  (1-6)  and  (1-9)  we  conclude  that  RZ  signaling  of  Type  II  is  3  dB 
superior  to  R  Z  signaling  of  Type  I. 

(b)  NRZ-L  or  Manchester  (Split-Phase)  Signaling.  NRZ-L  signals 
are  characterized  as  the  negatives  of  each  other,  in  fact,  * 


s  (t)  =  -s  (t)  O^tsT 


(I- 10) 


The  BEP  for  a  perfectly  synchronized  receiver  is  given  by 


P£  =  i  erfc(^) 


(1-11) 
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which  says  that  the  NRZ-L  signaling  format  is  3dB  superior  to  the 
Type  II  -  RZ  signaling  format.  This  BEP  can  also  be  obtained  with 
Manchester  coding. 

(c)  NRZ-M  and  NRZ-S  Signaling.  NRZ-M  and  NRZ-S  formats 
represent  systems  which  employ  differential  encoding  and  differential 
decoding.  For  these  ct-oes  the  BEP  can  be  shown  to  oe  given  by 

PE  =  2  exp(-R)  (1-12) 

For  R  >  10  dB,  NRZ-L  and  NRZ-M  or  NRZ-S  yield  approximately  equivalent 
bit  error  probabilities. 

(d)  Miller  Coding  (Delay  Modulation).  To  detect  the  Miller  code 

a  slightly  different  receiver  mechanization  of  the  data  detector  is  required. 

If  one  examines  only  the  data  received  during  a  single  bit  interval  then  the 
optimum  data  detector  is  illustrated  in  Fig.  1-3.  For  convenience  we 


qi  ~q2 


l  it; tire-  1-3.  Optimum  Sinqle-Sliot  Data  Detector  for  the  Miller  Code. 
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assume  that  the  transmitted  signals:  are  sequences  of  xl's.  For  this 
case,  the  BEP  is  easily  shown  to  be  given  by  (Ref.  3) 


P^,  =  erfc 


(1-13) 


For  R  »  1,  we  have  from  (1-13)  that 


„  ,  TP 

PE  *  erfc  <!  7 


d-14) 


so  that  the  BER  is  twice  as  large  as  that  obtained  with  Type  II  -  RZ 
signaling.  If  the  performance  of  Miller  codes  is  compared  with  NRZ-L 
'Manchester  coding  or  split  phase)  then  we  see  from  (1-11)  and  (1-14) 
that  the  BEP  for  NRZ-L  signaling  is  |  as  large  at  3  dB  less  signal-to- 
noise  ratio  as  the  BEP  for  Miller  coding. 

It  is  interesting  to  compare  the  performance  of  a  Miller  code  with 
quadriphase,  i.  e.  ,  four  phase  PSK  signaling.  For  quadriphase  signaling 
it  is  easy  to  show  that  (Ref.  3) 


P 


E 


erfc  «/R/2‘  ; 


7  eric 


R  »  1 


(1-15) 


so  that  for  R  »  1  the  BEP  for  quadriphase  signaling  is  equivalent  to 
the  performance  of  the  Miller  code. 

In  passing  we  note  that  it  appears  possible  to  derive  the  optimum  bit 
detectors  for  the  Miller  code  when  the  infinite  past  and  infinite  future 
of  the  data  are  considered  in  the  detection  process.  The  BEP  for  the 
rather  complicated  data  detectors  has  been  evaluated  and  it 
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has  been  shown  that  such  detectors  perform  superior  to  that  of  Fig.  1-3 
by  less  than  0.  5  dB.  The  question  as  to  whether  one  should  mechanize 
the  optimum  data  detector  for  Miller  coding  must  therefore  be  determined 
by  the  complexity  of  data  processing  allowed  in  at  the  point  of  reception. 
1-2.  BEP  Comparisons  for  the  Case  of  Perfect  Bit  Synchronization.  In 
summary  we  have  that  the  BEP  for  RZ,  NRZ  and  Miller  code  signaling  is 
given  by 

BERrz  r:  i  erfc(./R7?)  ;  Type  I-RZ 

EERrz  =  1  erfc(v/R7?)  ;  Type  II-RZ 

(1-17) 

be1Wz-l  * 

(Manches*e- 

BERMillrr  '  erfc(^?)[l  -  1  erfc(7R72)l 


Fig.  1-4  represents  a  plot  of  the  BEP,  as  determined  from  these  expressions, 

versus  the  signal-to-noise  ratio  R  =  ST/N^  =  E/N^  in  dB.  Thus,  using 

Miller  coding,  we  see  that  for  large  R  the  BEP  is  approximately  twice 

as  large  as  can  be  obtained  with  NRZ-L  with  one  half  the  power.  For 

-3  -5 

BEP's  on  the  order  of  10  or  10  Miller  coding  requires  approximately 
3.  5  dB  more  power  than  NRZ-L  or  Manchester  coding.  On  the  other 


hand.  Type  II-RZ  signaling  is  superior  to  Miller  coding  in  this  vicinity 
of  BEP  by  approximately  0.  5  dB  while  Type  I-RZ  is  inferior  to  Miller 
coding  by  2.5  dB. 

-3 

Finally,  for  R  ^  1  or  P  <10  we  note  that 

hj 


(1-18) 


where  the  symbol  4-PSK  implies  quadriphase. 


II.  SPECTRAL  OCCUPANCY,  BANDWIDTH  CONSIDERATIONS.  AND  THEIR  EFFECTS  ON 
BIT  SYNCHRONIZATION  (REFERENCE  3) 

Aside  from  BEP  considerations,  a  number  of  other  factors  must  be 
taken  into  consideration  when  attenr\pting  to  select  a  particular  baseband 
signaling  and  bit  synchronization  technique.  These  include: 

(1)  Signal  Bandwidth 

(2)  Signal  Bit  Synchronization  Capabilities 

(3)  Signal  Interference  and  Noise  Immunity  Properties 

(4)  Cost  and  Complexity  of  Implementation  of  the  Data  Detector 
and  Bit  Synchronizer 

In  what  follows  we  discuss  the  spectral  occupancy  of  the  baseband  signaling 
schemes  by  presenting  the  power  spectral  density  functions  and  discuss 
their  relationship  to  the  problem  of  bit  synchronization.  We  shall  define 


A  P 

S.(f)  =  s.(t)exp(-j2nft)dt 
1  J0  1 

to  the  Fourier  transform  of  the  elementary  signal  s.(t). 

(a)  Power  Spectral  Density  for  Type  I-RZ  Signaling.  Here  we  have 
a  situation  where  (f )  =  0  and 


-  n  TVS’  [" sin  TTfT/2*] 

Si'f)  =  2  If  TTfT/2 "]  ^PHTTfT/2) 


(II- 1) 


For  a  purely  random,  equal  probable,  source  the  power  spectrxim  is 
given  by 


S(f) 

E 


RZ 


6(f)  1 

'  16T  16T 


_1_ 

16 


n--°» 


/ sin  TTfT /2 
\  TTfT/2  / 

(II- 2) 
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Notice  the  delta  function,  6(f),  appearing  in  the  spectrum.  The  spike  at 
f  =  1/T  can  be  used  to  establish  bit  synchronization;  however,  the  signal 
energy  at  f  =  1/T  would  produce  a  self-jamming  effect  on  the  establishment 
of  bit  sync. 

(b)  Power  Spectral  Density  for  NRZ-L  Signaling.  Fo’*  NRZ  signals 
defined  by  (I- 10)  with  s^(t)  =  JS  for  0  £  t  £  T  then  the  power  spectral  density 
is  given  by 


Notice  that  there  exists  no  spikes  in  the  spectrum  and  that  S(0)  =  E.  These 
facts  are  of  considerable  importance  when  the  bit  synchronization  problem 
is  considered. 

(c)  Power  Spectral  Density  for  Manchester  Coding  (Split  Phase). 

In  this  case  it  is  easy  to  show  that 


S(f) 

E 


Manchester 


[sin(TTfT /2)1 4 
fnfT/2]  2 


(H-4) 


Notice  here  that  the  spectrum  contains  no  delta  functions;  however, 

S(0)  =  0.  In  contrast  to  RZ  signaling,  these  facts  make  Manchester  coding 
attractive  from  the  viewpoint  of  obtaining  bit  synchronization  from  the 
received  signal. 
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(d)  Power  Spectral  Density  for  Miller  Coding.  For  Miller  coding 


the  power  spectral  density  can  be  shown  to  be  given  by  (Ref.  1) 


S(f) 

E 


Miller 


2x  [17  +  8  cos  8x] 


[23  -  2  cos  x  -  22  cos  2x 


-  12  cos  3x  +  5  cos  4x  +  12  cos  5x 

+  2  cos  6x  -  8  cos  7x  +  2  cos  8xJ  (II- 5) 

where  x  =  rrfT.  Notice  here  that  S(0)  =  0.1  and  that  the  spectrum  contain 
no  spikes.  These  spectral  features  make  Miller  coding  attractive  from 
the  bit  synchronization  point  of  view. 

(e)  Spectral  Comparisons  for  RZ,  NRZ,  and  Miller  Coding.  Plotted 
in  Fig.  II- 1  are  the  power  spectral  densities  of  RZ,  NRZ,  Manchester 
coding  and  Miller  coding.  The  Spectral  properties  of  Miller  coding  which 
make  it  attractive  for  magnetic  tape  recording  (as  well  as  phase-shift 
keyed  signaling)  include: 

(1)  The  majority  of  the  signaling  energy  lies  in  frequencies  less 
than  one  half  the  bit  rate  %  -  1/T. 

(2)  The  spectrum  is  small  at  f  =  0.  This  spectral  minimum  facilitates 
the  problem  of  bit  synchronization. 

(3)  A  reduced  spectral  density  in  the  vicinity  of  f  =  0  is  important  in 
tape  recording  because  of  their  poor  d-c  response. 

(4)  As  a  result  of  (3)  a  lower  magnetic  tape  recording  speed  (higher 
packing  density  can  be  used). 

(5)  The  Miller  code  is  insensitive  to  the  180°  phase  ambiguity  that 
is  common  to  NRZ-L  and  Manchester  coding. 
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(6)  The  bandwidth  required  for  Miller  coding  is  approximately  | 
that  required  by  NRZ-L  and  ^  that  required  for  Manchester 
coding,  see  Fig.  II- 1. 


0  0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6  1.8  2.0 

fT 


Figure  II— I .  Spectral  Density  of  Random  NRZ, 
Manchester  Coding,  and  Miller  Coding. 


111.  PERFORMANCE  OF  PCM  SIGNAL  CONDITIONERS  IN  THE  PRESENCE  OF  BIT 
SYNC  JITTER  (REFERENCE  3) 

In  this  section  of  the  report  we  consider  the  combined  problem  of 
data  (bit)  detection  and  the  effects  of  bit  sync  jitter  on  BEP.  Our  method 
of  approach  is  to  evaluate  the  BEP  conditioned  upon  the  fact  that  the  bit 
synchronizer  under  consideration  produces  sync  pulses,  i.  e.  ,  samples 
the  output  of  the  data  detector,  which  are  in  error  by  t  second.  Such  a 
BEP  will  be  called  the  conditional  bit  error  probability  (CBEP),  We  shall 
denote  the  CBEP  by  P  (t).  This  approach  has  the  advantage  that  one  can 
consider  the  effects  of  bit  sync  jitter  on  BEP  for  a  broad  class  of  bit 
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synchronization  systems  without  having  to  specify  the  BEP  for  each 
situation.  The  average  BEP  is  then  obtained  by  averag  ng  the  CBEP  over 
the  probability  density  function  p(T)  of  the  bit  sync  jitter.  The  probability 
density  function  p(r)  is  determined  by  the  particular  bit  synchronization 
system  under  consideration.  Thus,  in  general,  we  have  that  the  average 
BEP  is  given  by 

PE  =  I  P<T)PE(T>dT  (m-1) 

T 

We  now  present  the  conditional  BEP  expressions  for  RZ,  NRZ,  Manchester 
and  Miller  coding. 

(a)  Conditional  BEP  for  Type  I  -  RZ  Signaling.  If  one  establishes  a 
threshold  at  ST/2  in  the  data  detector  and  assumes  that  the  bit  sync 
sampling  pulses  are  in  error  by  t  seconds  then  the  conditional  BEP  is 
easily  shown  to  be  given  by  (Ref.  3) 


pe(t) 


(HI- 2) 


This  result  assumes  rectangular  on-off  pulses  characterized  by  (I-  5  ). 
When  t  =  0  this  result  reduces  to  (I-  6  )  as  it  should. 

(b)  Conditional  BEP  for  Type  II  -  RZ  Signaling.  For  this  case  the 
signal  pulses  are  characterized  by  (I-  8  )  and  by  establishing  the  optimum 
data  detection  threshold  of  zero,  the  CBEP  is  easily  shown  to  be  given  by 


(Ref.  3) 


PE(M  =  4  \  erfc 


*  «"■=  =  hi  4 


(III  -  3 ) 
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when  the  bit  sync  error  is  t  seconds.  When  r  =  0  this  result  reduces 
to  (I-  9  )  as  it  should. 

(c)  Conditional  BEP  for  NRZ-L  S:?naljng.  The  signal  pulses  are 
characterized  by  (1-10).  For  a  t  second  bit  synchronization  error  the 
CBEP  can  be  shown  to  be  given  by  (Ref.  3) 


pe(t) 


\  <erfc(,/f?)  +  erfc 


(III-4) 


Notice  that  when  T  =  0  this  result  reduces  to  (I- 11  )  as  it  should. 

(d)  Conditional  BEP  for  Manchester  Coding.  Since  Manchester  coding 
requires  a  transition  during  every  bit  interval  the  CBEP  is  slightly  higher 
than  that  of  NRZ-L.  In  fact,  it  is  easy  to  show  that  (Ref.  3) 


pe(t) 


(HI- 5) 


when  the  bit  sync  pulse  is  in  error  by  T  seconds 

(e)  Conditional  BER  for  Miller  Coding.  Evaluation  of  the  CBEP 
for  Miller  coding  is  slightly  more  difficult  than  the  other  PCM  baseband 
signaling  techniques  considered  thus  far.  The  approach  to  developing 
the  theoretical  expression  for  the  CBEP  is  to  considei  the  data  stream 
as  composed  of  signal  vectors  s  =  (1,  1),  s  =  (-1,  1),  s  =  (-1,  -1)  and 

L  Ct  j 

s4(l,  -1).  The  conditional  BEP  given  that  s^  was  transmitted  is  given  by 

P  (t  j  s  )  =  P  (t  |  s  ,  no  transition)Pr[No  transition] 

XLi  I  Jii  I 

+  P  (T  |  s  ,  transition)Pr[transition]  =  Pt,(t  j  s  ) 

XL#  J.  -1  j 

(HI-6) 
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Now  the  joint  probability  density  function  of  the  correlator  outputs 
of  Fig.  1-3  given  that  Xj  is  received,  no  transition  occurs  and  the  bit 
sync  sampling  pulse  is  in  error  by  r  seconds  is  given  by 


l  r  (qj  -  >/m2  1  r  <q, 

p(q1,q2|s1,T.NTr)  =  —  exp^- - - - J  exp^-  — 


(UI-7) 


and  when  a  transition  occurs  we  have 


p(q1.q2l  sj» Tr)  = 


h  exp[- 


(qx  - 


(q2  -  JR  (l  -  ihl 


(III- 8) 


Using  these  p.  d.  f. 1  s  we  have 


PE(t!V  =  J  J  P(qi.q2!s1»T-NTr)dq1<iq2  +  f  [  p(q1,q2|s1,T,Tr)dq1dq2 

0  -®  ”0 

(m-9) 


This  can  be  reduced  to 


pe(tIV  -  pe(t|s3>  x  (/¥)[•  -  i  "*(/!)]  (III-10» 

t  i  erfc[/f  (l  -  -  I 


If  we  transmit  s^  =  (-1,1)  then  the  conditional  BEP  can  be  written  as 
P  (t|s  )  =  P  (t  I  s  )  =  P  (t|s  NTr)Pr(NTr|s  )  +  P  (T|s  Tr)Pr[Tr|s  ] 


2'  '  E'  1  V  ’  L  1  2- 
(III- 11) 
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If  we  note  that  Pr[Tr  ]  "1  =  0  and  that  P  (t  |  s  )  =  P  (t  |  s  .)  due  to  the  fact 

that  the  Miller  cods  is  insensitive  to  the  180  degrees  phase  shift  then  we 
can  write 

PE(T^S2)  =  PE(TK)  =  pEMs2’NTr>  (HI-12) 


This  probability  can  be  evaluated  in  terms  of  the  integral 

00  09 

PE(t!s2)  =  2J  1  P<Vq2,T’S2)dqidq2 


(III- 13) 


0  0 


where  the  joint  p.d.f.  pfq^q^lT,  s^)  is  characterized  by 


.  !  f  (qj+b)  1  f  (q2-b)  I 

P^1-<12!t*82)  =  2tt  eXP  2 - ^JeXP[' - 2 - 


(HI-14) 


Direct  substitution  and  simplification  leads  to 


pe(X)  =  I  +  1T  erf2(a)  +  ^  erfZ(/l)‘ 


4^  c  expf-fx-VR)^]  erf  xdx  - 


where 


3  pv  -  ~  2 

exp[  {-  (x-z/R^l  -4  |  \  |)}  ]  erf  xdx  (III- 15) 


a  =  Jf  0-4fXj)  ;  b  =  J~—  (1  -2 1 X | )  (HI-16) 


If  the  bit  sync  error  t  is  zero,  then  the  above  result  reduces  to  1-13 
as  it  should. 
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Ill- 1 -  Numerical  Evaluation  of  the  BEP 


In  order  to  numerically  evaluate  the  average  BEP,  it  is  convenient 
to  characterize  the  probability  density  function  p(T)  of  the  bit  sync  jitter 
as  a  function  of  the  normalized  bit  sync  error 

X  =  t/T  (HI-17) 

We  shall  postulate  a  probability  density  function  of  the  bit  sync  error  in 
such  a  manner  that  the  numerical  results  obtained  are  in  terms  of  the 
rms  normalized  bit  sync  jitter  a  .  This  approach  has  the  advantage  that 

A 

the  set  of  curves  obtained  of  somewhat  universal  in  that  the  effects  of  bit 
sync  jitter  on  BEP  is  a  function  of  the  rms  bit  sync  jitter,  say  Oy  is  clearly 
manifested.  Obviously 

aT=0XT  (III- 18) 

We  shall  assume  that  the  average  clock  frequency  of  the  bit  synchronizer  is 
equal  to  that  of  the  data  stream.  The  uncertainty  is  due  entirely  to  the  bit 
sync  jitter. 

In  terms  of  the  normalized  bit  sync  error  X=  t/T  the  average  BEP 
can  be  written  as 

P„  =  !  P_(X)p(X)dX  (III- 19) 

E  ,x  E 

where  p(X  )  is  the  probability  density  function  (p.  d.  f.  )  of  the  normalized 
error  X  =  t/T.  This  p. d.f.  is  conveniently  characterized,  in  terms  of 
,  by  (Ref.  2) 
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when  RZ  and  NRZ  signaling  is  employed.  Here  Iq(x)  is  the  modified 
Bessel  function  of  order  zero  and  imaginary  argument. 

For  PCM  baseband  signaling  techniques  which  employ  transitions  in 
the  middle  of  every  bit  interval  then  the  p.d.f.  of  the  normalized  bit  sync 
error  is  characterized  by 


Substitution  of  (III-2Q  or  (III-21  )  into  (III- 19)  with  the  appropriate  expression 
for  P  (X)  and  performing  numerical  integration  on  a  digital  computer 

hj 

leads  to  the  graphical  results  of  Figs.  HI-1,  HI-2,  HI-3  and  III-4  for  BER 
vs.  the  signal-to-noise  ratio  R  =  ST/N^  with  the  normalized  bit  sync 
rms  jitter  0^-  0_/T  as  a  parameter. 
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IV,  BEP  FOR  NRZ-L  SIGNALING  WITH  ARBITRARY  WAVEFORMS 


We  now  generalize  the  previous  theory  to  include  conditions  where 

Gaussian  poise  is  added  to  pulses  of  arbitrary  shape.  Such  situations 

arise  when  the  PCM  data  stream  is  filtered  prior  to  detection.  As  before 

we  consider  the  situation  where  the  bit  sync  is  in  error  by  T  =  XT  sec. 

th 

The  cross  correlator  output  on  which  a  decision  is  made  for  the  n 
transmitted  bit  now  depends  on  the  received  data  in  the  time  interval 

A  A  A 

(n-l)T  l  e^t  s  nT  +  e  where  e  is  the  bit  timing  estimate  obtained  for 
the  bit  synchronizer.  Stated  mathematically. 


where  a^  is  the  n  detected  bit  and  pg(t)  is  the  basic  pulse  shape  of  the 
PCM  signal.  The  received  PCM  data  stream  plus  noise  is  characterized 
by 

K 

x(t)  =  ^  akpg[t  -  (k-l)T  -  e]  +  n(t)  (IV-2) 

k=0 

th 

and  a^  =  ±1  during  the  k  transmission  interval.  Here  e  serves  as  the 
pulse  epoch  or  initial  phase  and  K  represents  the  number  of  bit  intervals 
that  e  is  considered  to  be  constant.  If  this  expression  is  substituted  into 
the  above  expression  for  a^  then  by  breaking  Q  into  the  signal  component 
Qg  and  noise  component  Q  then  it  is  easy  to  show  that  (Ref.  3) 

Qs  =  an-lRsC'T(1-X)1+anRs(XTi  +  an+lRs^T(1  +  X)]  ^-3) 


[t  -  (n-l)T  -  € ]dt  /  =  sgn  Q 


(IV-l) 


=  sgn 


!r 


nT+e 


iJ  .  X(t)ps 

Un-l)T+e 
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where  R  (t)  is  the  signal  pulse  autocorrelation  function  R  (t)  defined  by 

8  S 

.  rT-T 

R  (T)  =  ]  p  (t*+T )p  (t')dt'  (IV -4) 

s  8  8 

With  regard  to  the  input  bits  affecting  the  calculation  of  Q  ,  one  of  four 

s 

th 

situations  can  occur:  either  the  n  bit  and  its  adjacent  bit  are  alike  or 

they  are  unalike  (negative  followed  by  positive,  or  vice  versa).  In  fact,  if 

X  >  0  then  a  and  a  ,  effect  the  decision  while  if  X  <  0  then  a  and  a  , 
n  n+1  n  n-1 

enter  into  the  decision  variable.  Hence,  the  CBEP  is  given  by  (Ref.  3) 

2 

Pg(X)  =  ^  ^  [Prfn**1  and  adjacent  bits  =  (-l)k}FV{error  in**1  bit  =  (-l)k} 
k=l 

+  Pr{nth  bit  =  (-l)k  and  adjacent  bit  =  {-l)k+1]  *  (IV-5) 

Pr {error  |n4^  bit  =  (-l)k}] 


Since  the  random  variable  Q  is  conuitionally  Gaussian  on  a  ,  ,a  ,  a 

’  n-1  n  n+1 

and  t  =  XT  fixed,  the  input  bits  occur  with  equal  probability  and  are 
2 

independent,  and  Oq  =  N^E/2,  it  is  easy  to  show  that  (Ref.  3) 

P~(t)  =  jerfck/R^r  (t)  +  r  (1  -  |t|)]} 

S  S 

+  i  erfcWR[r  <t)  -  r  (1  -  |t  |)]}  (IV-6) 

8  S 


where 


/  X  A 
rs(T)  = 


Rs(T) 

E 


The  average  BEP  is  found  from  (III- 19)  and  (III- 20 >  using  this  result. 


27 


In  passing  we  note  that  for  rectangular  pulses 


r  _(t) 
s 


M  s  T 
M  *  T 


(IV  -7) 


To  obtain  numerical  results  for  various  filtered  PCM  signals  numerical 
integration  is  required. 


V.  BIT  SYNCHRONIZATION  SYSTEM  PERFORMANCE  MEASURES  (REFERENCE  2) 

In  this  section  of  this  report  we  shall  digress  from  the  problem  of 
bit  detection  and  consider  the  compenion  problem  of  bit  synchronization. 

The  motivation  here  is  very  simple:  There  has  been  very  little  useable 
theory  developed  which  allows  one  to  determine  the  tradeoffs  which  can 
be  realized  in  the  design  and  operation  of  a  bit  synchronization  system 
(BSS)  during  bit  sync  acquisition  and  bit  sync  ti  eking.  This  is  particularly 
true  when  the  input  to  the  BSS  is  a  PCM  signal.  The  lack  of  concentration 
on  formulating  and  characterizing  the  fundamental  limitations  of  BSS 
has  lead  to  difficulties  in  not  only  procuring  BSS  but  also  in  implementing 
them  such  that  accurate  and  predictable  procurement  specifications  can  be 
met.  For  example,  builders  of  commercial  bit  synchronizers  usually  do 
not  adequately  specify  or  characterize  the  acquisition  or  tracking  perform¬ 
ance  of  a  particular  bit  synchronizer.  Certainly,  the  state  of  affairs  is 
such  that  tradeoffs  in  the  two  modes  (tracking  and  acquisition)  cannot  be 
satisfactorily  made.  For  example,  a  BSS  configuration  is  usually  changed 
by  varying  the  system  "bandwidth"  by  means  of  a  switch  on  the  front  panel. 
The  primary  reason  for  this  situation  is  that  the  theory  has  not  been 
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satisfactorily  advanced  oweing  to  the  nonlinear  nature  of  a  multi¬ 
dimensional  problem. 

One  purpose  of  this  ,'eport  is  to  discuss  the  many  parameters  which 
serve  to  characterize  the  performance  of  a  BSS.  This  discussion  is  given 
so  that  people,  who  are  responsible  for  procuring  BSS,  can  become 
aware  of  what  is  required  in  writing  a  bit  sync  specification,  what  parameters 
are  reasonable  and  meaningful  to  spec  and  what  constitutes  physically 
realizable  specifications.  Before  attempting  to  set  forth  any  theory  it 
seems  appropriate  to  discuss  the  many  facets  of  the  bit  sync  problem. 

This  will  set  the  tone  for  that  which  follows. 

V  -1.  Operational  Behavior  and  Performance  Characterization  of  bit 
Synchronizer  s.  The  operation  of  any  bit  synchronization  system  (BSS) 
is  conveniently  discussed  by  considering  its  two  modes  of  operation, 
namely,  the  acquisition  mode  (achieves  the  synchronous  state)  and  the 
synchronous  or  tracking  mode  (retains  the  synchronous  state).  On  the 
other  hand,  the  performance  of  any  BSS  is  characterized  by  certain 
performance  measures  which  are  statistical  in  nature  when  noise  is  present 
and  deterministic  when  noise  is  absent.  In  order  to  introduce  these 
performance  measures  we  shall  discuss  the  operational  behavior  of  BSS 
from  the  viewpoint  of  rotating  phasors.  This  point  of  view  has  the  distinct 
advantage  of  being  simple  without  losing  generality  in  illustrating  the  oper¬ 
ational  behavior  of  a  BSS. 

A 

Let  the  real  part  of  the  phasors  exp[i$(t)"|  and  exp[i$(t)"j  be  associated 

A 

with  the  signals  s(t,  $)  and  r(t,  $),  respectively.  These  phasors,  which 

A  ' 

rotate  with  instantaneous  angular  velocities  uu(t)  =  $(t)  rad/sec  and 
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u?(t )  =  $*' • ,  j  ospectively,  are  shown  in  Fig-  V-l.  The  bit  sync  error, 

A  A 

cp(t)  =  ${ti  -  ?(t),  has  temporal  characteristics  which  are  strongly  dependent 
on  the  Application  and  mode  of  operation;  however,  the  performance  measures 
associated  with  the  random  process  fPft)}  are  essentially  independent  of 
the  application.  We  note  that  when  the  phase  function  $(t)  is  characterized 
by  a  random  process  the  signal  phasor  expfi$(t)"]  and  reference  phasor 

A  ^  A 

exp[i$(t)]  rotate  randomly  at  rates  proportional  to  $ (t)  and  $(t). 


I'  il’iire  \-l.  I’lmsor  hiui’rani  lor 
Describing  Operation  of  a  HSS. 


(a)  Acquisition  Mode.  Suppose  the  system  enters  the  acquisition  mode 

_  A 

at  t  -  t^.  Initially  the  phasors  exp[i$(t)"|  and  exp[i$(t)]  rotate  at  angular 
velocities  &  =  'JJftg)  and  =  ift^)  rad/ sec  respectively.  Thus,  initially 
exp[il(t)1  rotates  at  rate  =cp(t^)  =  =  ou  -  rad/sec  relative  to  exp[i£(t)"l. 

The  parameter  Q^/2tt  is  called  the  initial  frequency  detuning  in  the  system. 

Figure  V-2  shows  a  typical  behavior  of  cp(t)  and  cp(t)  for  a  second-order 
BSC  operating  in  the  acquisition  mode  without  an  acquisition  aid  in  the 
absence  of  noise.  From  this  figure  it  is  clear  that  the  behavior  of  Cp(t) 
and  cp(t)  during  acquisition  is  highly  dependent  upon  the  initial  state 
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Figure  V-2.  Typical  Behavior  of  <j>  and  <£  for  a 
BSS  in  the  Acquisition  Mode  (Noise  Absent). 


Xq  =  (t q) ,cp (^q)  1  =  fp0,CP0l*  ^et  synchronous  mode  be  defined  by  the 
conditions  |cp { t )  |  £  e  ,  |cp(t ) -2nTT |  ,  n  any  integer,  for  t  =  t  .  Suppose 

3. 

these  conditions  are  satisfied  at  t  =  t  -t  for  the  first  time.  The 

acq  a  0 

parameter  t&  is  called  the  signal  acquisition  time  (often  called  pull-in 

time  or  lock-in  time)  since  it  is  the  time  required  to  reach  the  synchronous 

mode.  When  noise  is  present  the  situation  is  much  more  complicated  and 

we  do  not  discuss  all  the  details  here;  however,  for  this  case  t  become 

acq 

a  random  variable.  Thus  an  important  performance  measure  for  the 

acquisition  mode  is  the  statistical  moments  Tn  =  E(tn  );  E'-  )  denotes 

acq  acq 

expectation  in  the  statistical  sense.  Obviously,  depends  upon  the 

initial  state  of  the  system.  When  the  initial  state  is  random  and 

characterized  by  a  probability  density  function  (p.d.f, ),  then  T°  must 

acq 

be  averaged  over  this  density. 
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Shortly  after  acquisition,  say  t>t^  seconds,  the  "steady- state"  mean 
<P/2tt  =  (ti)  -  o)^.)/2tt  represents  the  mean  residual  frequency  detuning.  Here 
w  is  the  average  radian  frequency  of  the  synchronized  generator  and 
'JUy)/2TT  is  the  mean  frequency  shift  of  the  waveform  generator  to  be 
synchronized.  When  l(t)  =  <t'  is  constant,  the  largest  value  of  say 
|fiQ|  ,  from  which  the  system  can  reach  the  synchronous  mode  is  called 
the  signal  acquisition  (often  called  pull-in  or  lock-in)  range;  that  is,  l^glm 
represents  the  maximum  relative  angular  velocity  between  phasors  exp[i$(t)] 

A 

and  exp[i$(t)]  for  which  acquisition  is  possible.  Since  noise  is  present, 
characterization  of  the  probability  of  acquisition  at  time  t  is  also  of  interest. 

In  the  acquisition  mode,  a  BSS  is  highly  nonlinear.  In  fact,  the  so- 
called  linear  theory  cannot  be  used  to  carry  out  a  particular  design  nor 
can  it  be  used  to  account  for  the  system's  operational  behavior;  that  is, 
the  acquisition  problem  is  characteristically  nonlinear. 

(b)  Synchronous  or  Tracking  Mode.  Once  the  system  is  in  the 
synchronous  mode  a  more  or  less  steady- state  condition  exists  and  the 
system  tends  to  automatically  make  adjustments  so  as  to  maintain 
|cp(t)-2rinl^  e  and  jcp(t)  |  with  high  probability;  that  is,  it  tends  to 
maintain  the  synchronous  state.  Due  to  the  effects  of  noise,  n.(t),  and 
instabilities  in  the  generator  to  be  synchronized,  random  fluctuations 
in  cp(t)  take  place;  see  3. 
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In  the  synchronous  mode  a  number  of  performance  measures  are  required  to 

characterize  and  explain  system  behavior.  The  reason  is  that  the  performance 

measures  which  characterize  this  mode  of  operation  are  complicated  to  account 

for  when  ot.h  the  system  nonlinearity  and  noise  are  taken  into  consideration. 

If  the  system  is  linearized  to  simplify  the  problem,  the  performance  measures 

which  sei  ->  to  characterize  the  fundamental  behavior  usually  have  no  meaning; 

that  is, operation  in  the  synchronous  mode  is  also  characteristically  nonlinear. 

In  this  mode  the  time -dependent,  conditional  transition  p.d.f. 

P(cp(t!Vto)  of  the  absolute  phase  error  cp(t)  is  of  interest.  Furthermore, 

the  reduced  phase  error  ^^(t)  =  [cp(t)  +  ttlmod  2TT  +  (2n-l)TT  (such  as  displayed 

by  a  phase  meter)  is  of  interest.  Here  n  is  considered  to  be  any  fixed 

integer.  See  Fig.  V-4  for  a  typical  sample  function  ofcp^(t)  and  cp(t).  )  In 

loose  terms,  the  bit  sync  error  process  [cp  (t)}  is  called  the  modulo-2TT 

r 

process.  Thus  the  time-dependent,  conditional  transition  p.  d.  f.  p(cp  ,  t  t^,  n) 
of  the  bit  sync  error,  cp ^(t ),  which  a  phase  meter  would  read  and  the 
steady-state  conditional  transition  p.d.f.  pft>  |n)  are  of  importance. 

With  the  passage  of  time,  noise  causes  bit  slippage  to  occur.  Here 
we  define  bit  slippage  as  the  change  of  Cp  from  cp^  tocp^  +  2tt  ortpQ  -  2tt,  i.  e. , 
ACP/2TT  ±1  cycle;  see  Fig.  V-4a.  This  corresponds  to  the  situation  in 
Fig. V-4  where,  for  any  initial  state  =  [^q'^Pq]*  the  phasor  exp[i$(t)J 
rotates  ±2tt  radians  relative  to  the  phasor  exp[i$(t)].  In  physical  terms 
the  generator  to  be  synchronized  drops  or  adds  one  cycle  of  oscillation 

A 

in  r(t,§)  relative  to  s(t,  §);  or  stated  another  way,  the  signal  s(t,  $)  is 

A 

is  shifted  ±T  seconds  relative  to  the  reference  signal  r(t,$).  This 
phenomenon  is  called  bit  slippage. 
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its  p.  d.  f.  or  its  moments,  i.  e.  ,  E(T  ).  The  first  moment  of  this  time 

fp 

is  called  the  average  time  to  first  loss  of  bit  sync  and  denoted  by  E(T^). 

'Willinm  <  Lindsey , Synchronization  Systems  in  Communication  ami  Control,  I c  1  1972.  Reprinted  by  permission  of 
Prentice-Hall,  Inc  ,  Knglewood  Cliffs,  New  Jersey 
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The  probability  of  bit  sync  failure,  P(t),  in  the  time  interval  [t  ,  t] 
as  well  as  the  probability  of  slipping  n  bits,  P(N-n),  in  the  time  interval 
[tg.t]  are  of  interest.  Since  biases  are  frequently  present  in  BSS,  the 
average  number  of  clockwise  (counterclockwise)  rotations  per  second  of 

A 

the  phasor  exp[i$(t)]  through  2rr  relative  to  the  phasor  exp[i$(t)]  and  the 
respective  probabilities  are  of  interest.  We  shall  refer  to  the  dropping 

A 

and  adding  of  oscillations  in  rft,  $)  relative  to  s(t,  $)  as  bit  slippage  to  the 
right  and  left,  respectively.  The  average  number  of  bits  slipped  to  the 
right  (left)  per  second  is  characterized  by  N^(N  ).  Other  parameters 
which  require  a  statistical  characterization  by  means  of  the  nonlinear 
theory  include  the  average  time,  AT,  between  bit  slipping  events  and  the 
average  time  the  BSS  remains  out  of  the  synchronous  state;  see  Fig.  V-4a. 

Table  V-l  summarizes  the  measures  of  performance  which  serve  to 
characterize  the  fundamental  behavior  of  BSS  systems.  It  is  clear  that  the 
design  or  procurement  engineer  has  a  multitude  to  consider;  however, 
in  what  follows  we  shall  be  primarily  concerned  with  bit  slippage  rate, 
bit  sync  jitter,  bit  sync  acquisition  time  and  bit  sync  acquisition  range. 

The  bit  error  probability  performance  has  been  considered  in  an  earlier 
part  of  this  report.  A  rather  complete  bibliography  of  reference  material 
is  given  in  the  list  of  references. 


Table  V-l.  Performance  Characterization  of  BSS  (Reference  2)* 


Acquisition  Mode 

Tracking  Mode 

1. 

Bit  sync  acquisition  range 

1. 

Steady  state  p.d.f.  of  bit  sync 

2. 

Bit  sync  acquisition  time 

2. 

Statistical  moments  of  the  bit  sync 

3. 

Probability  of  acquisition  at 
time  t 

3. 

Moments  cf  time  to  first  loss  of 
synchronization  (sync) 

4. 

Bit  sync  acquisition  figure  of 
merit 

4. 

Average  time  to  first  loss  of 
bit  sync 

5. 

Bit  sync  acquisition  behavior 

5. 

Average  number  of  bits  slipped 
per  unit  time 

6. 

Probability  of  synchronization  failures 
as  a  function  of  time 

7. 

Average  number  of  bits  slipped 
to  the  right  and  to  the  left  per 
unit  time 

8. 

Probability  of  k  slips  in  T  seconds 

9. 

Average  time  between  bit -slipping 
events 

10. 

Average  time  duration  of  failure 
events 

11. 

Mean  and  variance  of  the  reference 
generator's  frequency  error 

12. 

Phase  (frequency)  error  diffusion 
constant 

13. 

Signal  hold-in  range 

VL.  MODEL  OF  A  CLASS  OF  PCM  BIT  SYNCHRONIZATION  SYSTEMS 

Bit  synchronization  has  to  do  with  determining  the  instants  in 
time  when  the  modulation  may  change  states.  One  approach  to  pro¬ 
viding  this  timing  information  is  to  utilize  a  separate  communication 
channel  solely  for  establishing  bit  sync.  More  efficiently,  however, 
are  the  techniques  whereby  this  level  of  timing  information  can  be  ob¬ 
tained  directly  from  the  received  PCM  data  bearing  signal.  For  this 

*  William  C  Lindsey.  Synchronization  Systems  w  Cimimunicathm  and  Cnnimf,  <c)  1972.  Reprinted  by  permission  of 
PrenMce-llall.  Inc..  I  nplcwood  ClitK  New  Jersey 
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case  a  broad  class  of  bit  synchronization  systems,  suggested  by  maxi* 
mum  a  posteriori  estimation  techniques,  can  be  represented  by  the 
diagram  of  Figure  VI.  Specific  phase  detector  characteristics  in¬ 
cluded  are  the  early- late  gate  type,  those  which  incorporate  am  absolute 
value  approach,  a  difference  of  squares  approach  or  a  hybrid  of  these 
approaches,  amd  the  decision-directed  type.  We  shall  discuss  these 
approaches  later  on  in  this  section;  however,  before  we  get  involved 
with  detailed  comparisons  of  BSS  circuits  it  is  appropriate  to  discuss 
the  problem  on  the  basis  of  presenting  various  theorteical  results  which 
are  available  in  rather  general  terms. 


INPUT  - 
SIGNAL 


Figure  VI.  Model  of  a  Class  of  I5it  Svnclironizalion  Svsleins  (liSS). 


The  math  model  for  this  class  of  BSS  is  illustrated  in  figure  VI-1. 
Here  T(t)  represent  the  time  varying  epoch  to  be  tracked  while  x(t) 
represent  the  bit  sync  estimate  and  cp(t)  =  T(t)  -  r(t)  represents  the 
bit  sync  error.  The  phase  detector  characteristic  g(c p,  t)  is  usually 
time- varying  during  the  acquisition  mode  and  can  be  considered  to  ap¬ 
proach  the  steady  state  value  g(cp)  when  bit  sync  is  achieved.  The  con¬ 
stant  A  represents  the  rms  amplitude  of  the  bit  stream  K  is  a  gain 
determined  by  implementation  and  n  (t)  is  the  equivalent  additive  noise 
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process  which,  in  general,  depends  upon  g(cp,  t).  In  fact,  the  stochastic, 
intergro-differential  equation  of  operation  is  given  by 


cp  =  T (t)  -  — ^[Ag(cp,  t)  +  n  (t)] 
P  O 


where  p  =  d/dt  denotes  the  linear  operator  which  implies  differentiation 
with  respect  to  time.  In  fact,  the  above  equation  is  a  shorthand  way  of  writing 


dtp 

dt 


dt(t) 

dt 


AK 


-t 

f(t-X)[g(cp,  X)  +  n  (X)]  dX 

-00  ® 


V“ 


l-'igurt-  \  l-|.  Math  M(kIcI  for  a  Broad  Class  of  liSS. 


VI- I,  Acquisition  Time  and  Acquisition  Range  For  Two  Different  Phase 
Detector  Characteristics 


For  the  sake  of  definiteness  we  consider  BSSs  in  which  the  loop 
filter  assumes  the  form 


1  +  V  1  -  F 

F(s)  =fTTT  =  fo+TT7T 


r 


(VI-l) 


r 


where  Fq  =  T 2  / T j-  Such  a  filter  has  the  advantage  of  being  capable  of 
injecting  a  "quick  reaction"  type  component  of  voltage  directly  into  the 


38 


local  clock  generator  to  be  synchronized  as  well  as  possessing  (via  a 
storing  charge  on  a  capacitor)  memory.  For  this  case  a  class  of  second  - 
order  BSS  which  employ  the  phase  -lock  principle  are  represented  (in  no 
noise)  by  the  equation  of  operation 

3F=  °0  -  AKF0[g(y,]  *  y, 

(VI-2 

1  _  v 

'ris-'rrr^>tAKewJ 

Here  &f  =  0q/2tt  represents  a  frequency  offset  between  the  incoming  and 
local  clocks,  g(x)  represents  the  phase  detector  characteristic,  AK 
represents  the  bit  sync  loop  gain  and  x  represents  the  instantaneous  bit 
sync  error.  From  this  discussion  we  see  that  the  acquisition  time  is 
a  function  of  many  variables  including,  in  particular,  the  shape  of  the 
phase  detector  characteristic.  For  example,  as  we  set 

g(x)  =  sgn  [sin  x]  (VI-3) 


the  i  it  can  be  shown  (Ref.  2)  that  the  time  Tf  required  to  eliminate  bit 

X 

slippage  is  approximated  by 


(VI-4) 


where  Ai  0^/2^,  B^  is  the  single  side  loop  bandwidth,  and  r  =  4C 
where  Q  s  the  loop  damping  coefficient. 
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It  is  interesting  to  consider  the  case  where  the  phase  detector 
characteristic  is  sinusoidal,  i.  e. ,  g(x)  =  sin  x.  For  this  case  it  has 
been  shown  that  (Ref.  3) 


Comparing  this  result  with  that  above,  we  see  that  the  acquisition  time 
for  a  sync  system  in  which  g(x)  =  sgn  [sin  x]  is  approximately  3dP  faster 
than  that  of  a  sync  system  in  which  the  phase  detection  characteristic 
is  sinusoidal.  Fig.  VI-2  represents  a  plot  of  (VI-5)  that  is,  the  frequency 
acquisition  time  versus  | A f | /B,  for  a  Iood  damping  of  0.707  and  various 
values  of  F  . 

The  acquisition  range  | A f lrnax  is  also  interesting  to  compare  as  a 
function  of  the  phase  detector  characteristic.  In  fact,  for  g(x)  =  sgn  [sin  x], 
it  has  been  shown  that 


'F o(1  -  F0> 


(Vi-6) 


while  for  a  sinusoidal  phase  detector  characteristic 


AK 

Ztt 


(VI-7) 


Therefore,  when  Fq«  1,  the  acquisition  range  for  a  rectangular  phase 
detector  characteristic  is  Jz  times  larger  than  for  the  sinusoidal  phase 
detector  characteristic. 

Finally,  it  is  interesting  to  compare  the  time  to  slip  a  cycle,  say 
T2tt>  near  the  boundary  of  synchronization  for  various  phase  detector 
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characteristics.  Since  the  derivation  of  such  results  are  quite  lengthy 
we  merely  present  the  results  which  were  derived  in  (Ref.  2).  There 
it  is  shown  that  the  time  to  slip  a  cycle 

;  g(x)  =  sgn  [sin  x] 

;  g(x)  =  sin  x  (VI-8) 

;  g(x)  =  ^  |cp|  <  7T 

VI-2.  Bit  Slippage  Rate  as  a  Function  of  the  Phase  Detector  Characteristic 
The  bit  slippage  rate,  say  S,  in  a  PCM  bit  synchronizer  is  an 
important  parameter  in  that  it  gives  a  measure  of  the  rate  at  which  the 
bit  sync  reference  fails.  Although  this  parameter  is  difficult  to  account 
for  in  theory,  it  is  instructive  to  present  the  theory  to  date  for  those 
cases  for  which  it  is  analytically  possible  to  evaluate. 

If  one  considers  the  problem  of  obtaining  the  average  number  of 
bits  slipped  per  second  in  a  first-order  BSS  in  which  there  exists  zero 
frequency  detuning  and  white  noise  present,  then  it  is  possible  to  show 
that  the  normalized  slippage  rate  is  given  by  (Ref.  2) 


S_  _  _ (4/q) 


-  n  cp 

J  exp[  a  J  g(x)dx]  dcp 

-  TT  CP 

J  exp  [-  a  J  g(x)  dxjdcp 

-TT 

-  - 

-TT 

(VI-9) 


2 

where  a  =  [l/2nk0^  and  g(x)  represents  the  phase  detector  characteristic. 
Here  (X.  =  0  / T  represents  the  hit  sync  jitter  and  k  =  1  for  PCM  signaling 
formats  wherein  no  switching  takes  place  in  the  middle  of  the  bit  interval 
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and  k  =  2  for  PCM  signaling  techniques  in  which  switching  is  possible  in 
the  middle  of  a  bit  interval. 

For  the  case  where  there  exists  frequency  offsets  between  the  two 
clocks  then  it  can  be  shown  (Ref.  2)  that 


_s_  = 

WL 


4C, 


exp 


(VI-10) 


where  y  =  2nAf/AK  and  Cq  is  a  constant  which  normalizes  the  steady  state 
probability  density  function  of  tin-  bit  sync  jitter. 

Let  us  now  consider  some  special  cases  of  the  slip  rate  for  two 
special  phase  detector  characteristics.  Consider  first  the  phase 
detector  represented  by  g(x)  =  sgn  [sin  x].  For  this  phase  detector 
characteristic  it  is  easy  to  show  that  (Ref.  2) 


S  _  4  exp  (nq)  -  1 

W  ~  tt  na 

■L*  L 


(VI-11) 


while  for  the  sinusoidal  phase  detector  characteristic 


_S_ 

w. 


1 


n2al^(a) 


(VI-12) 


Here  W  =  2B  represents  the  double  sided  loop  bandwidth.  Figures  VI-3 
and  VI-4  represent  the  plot  of  S/2W^  for  a  sinusoidal  phase  detector 
characteristic  and  compares  this  with  the  S/W  obtained  for  the  rectangular 

I -J 

phase  detector  characteristic  g(x)  =  sgn  [sin  xj.  As  before  CT^  =  0  / T 
represents  the  normalized  bit  sync  jitter. 


43 


VII.  MATH  MODELS  AND  1 HEIR  RELATIONSHIP  TO  COMPUTER  SOFTWARE 

In  order  to  develop  a  computer  software  program  for  a  BSS  one 
must  determine  several  things.  First  of  all,  it  is  desirable  to  have  a 
baseline  BSS  which  is  optimum.  This  is  important  because  the  per¬ 
formance  measures  which  characterize  its  behavior  serve  as  a  limit 
beyond  which  no  other  BSS  can  outperform.  The  problem  of  selecting 
an  optimum  BSS  for  purposes  of  implementation  or  for  carrying  out 
performance  comparisons  largely  remains  an  open  question.  For 
example,  it  is  possible  to  use  the  techniques  of  maximum  a  posteriori 
(MAP)  estimation  and  maximum  likelihood  estimation  in  order  to 
specify  an  optimum  bit  sync  circuit.  Usually,  owing  to  the  nature  of 
the  statement  of  the  problem,  such  a  synchronizer  turns  out  to  be  an  open 
loop  device.  Furthermore,  the  criterion  of  selecting  the  peak  of  the 
a  posterior  probability  density  function  or  the  likelihood  function  does 
in  ro  way  stress  the  importance  of  the  acquisition  problem.  It  should 
also  be  noted  that  the  requirement  of  tracking  bit  sync,  once  it  is  obtained, 
places  even  stronger  requirements  on  the  MAP  or  maximum  likelihood 
estimate.  To  force  this  requirement  on  the  MAP  estimator  frequently 
leads  one  to  a  solution  in  which  the  MAP  estimator  is  physically 
unrealizable;  however,  the  open  loop  estimate  can  be  used  in  a  feedback 
loop  to  form  a  control  loop  which  will  track  bit  sync;  however,  the 
resulting  close  loop  synchronizer  which  is  suggested  by  the  MAP  estimator 
is  suboptimum.  These  approaches  have  been  the  subject  of  much 
analytical  and  experimental  investigation  (see  list  of  References). 

For  a  maximum  likelihood  synchronizer  Wintz  and  Luecke 
(Ref.  5)  evaluate  the  mean  magnitude  bit  sync  error  X  as  a  function  of 
the  PCM  signal  shape  and  a  parameter,  say  T  ,  which  serves  as  an 
estimate  of  the  number  of  bit  durations  that  bit  sync  remains  constant. 
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Quantitatively,  the  mean  magnitude  bit  sync  error  varies  inversely  as 


the  square  root  of  Tm  and  R  =  ST/Nq,  i.  e. , 

E I X I  T  R 
m 

The  results  of  Monte  Carlo  computer  simulation  of  performance  are 
illustrated  in  Figs.  VII-1  and  VII-2. 


Figwe  VII— I .  Mean  Magnitude  Bit  Syne  Krror  Performance  of 
the  MAI’  Synchronizer  as  a  Function  of  Memory  Time  (Raised 
Cosine  Pulses  Assumed).  (References  12.  13,  IS,  IS.) 

‘Reprinted  by  permission  of  the  Institute  of  Electrical  and  Electronics  Engineers.  Inc.,  345  East  47th  Street,  New  York,  N.Y. 
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Figure  V 11-2.  Mean  Magnitude  Bit  Syne  Krror  Performance  of 
the  MAP  Synchronizer  for  Three  Different  Pulse  Shapes  and  a 
Memory  Time  of  Fight  Symbol  Periods.  (Reference  5.) 


Other  BSS  suggested  by  closing  the  loop  using  a  MAP  estimator 
are  suggested  in  Figs.  VII- 3  to  5.  On  the  other  hand,  the  digital  data 
transition  tracking  loop  (DITL)  of  Fig.  VII. -6  represents  a  BSS  which 
employs  decision  feedback  and  an  error  channel  in  which  a  "window" 
can  be  used  to  improve  tracking  performance.  The  parameter 
is  used  to  adjust  the  window  width.  An  early-late  gate  symbol 
synchronizer,  its  associated  hit  detector,  and  three  different  phase 
detector  topologies  are  illustrated  in  Fig.  VII-7  and  8.  Finally,  an 
absolute  value  type  of  BSS,  say  AVTS,  is  illustrated  in  Fig.  VII-9. 

‘Reprinted  by  permission  of  the  Institute  of  Electrical  and  Electronics  Engineers,  Inc.,  345  East  47th  Street,  New  York,  N.Y. 
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Ki  i'll  re  V  II— S.  An  Karl> -I  -ale  Late  Type  of  Hit  Synchronizer  With 
a  Square  I -aw  Type  Vonlinearity  (Reference  .?)■ 


Figure  VII— 6.  The  Digital  Data  Transition  Tracking  Loop  Type 
Hit  Synchronizer  (l)TTL).  (References  12,  13,  IS,  IS.)* 

•Reprinted  by  permission  of  the  Institute  of  Electrical  and  Electronics  Engineers,  Inc.,  345  East  47th  Street,  New  York,  N.Y. 
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Kiijiire  MI-7.  Karly-I.atc  (kite  Hit  Synchronizer  and 
Associated  Bit  Detector.  (References  3,  IN.)* 


Reprinted  by  permission  of  the  Institute  of  Electrical  and  Electronics  Engineers.  Inc.,  345  East 
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Figure  Ml-S.  Three  Different  Phase  Detector  Topologies. 

(Hcferenees  3,  IS.)* 

•Reprinted  by  permission  of  the  Institute  of  Electrical  and  Electronics  Engineers,  Inc.,  345  Fast  47th  Street,  New  York,  N  Y. 
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VII- 1.  Performance  Comparison  of  Several  BSS  on  the  Basis  of  Bit  Sync 
Jitter  and  Slip  Time 


In  comparing  the  performance  of  several  different  BSS,  one  must 
choose  a  fixed  operating  condition  that  is  common  to  all.  One  possible 
basis  of  comparison  is  on  the  basis  of  equal  loop  bandwidths  for  all 
configurations  at  every  R  =  ST/N^.  This  comparison  is  meaningful 
when  considering  the  actual  design  of  a  loop  having  a  specified  band¬ 
width  operating  over  a  given  range  of  input  signal-to-noise  ratios. 

Fig.  VH-10  plots  the  variance  of  the  normalized  bit  sync  jitter 
0^  versus  R  with  6=  1/TW^  as  a  parameter.  Included  in  Fig.  VII-10 
are  the  corresponding  curves  for  the  DTTL  with  full  mid -phase  inte¬ 
gration  in  the  error  channel  and  the  DSL.  The  asymptotic  behaviors  of 
2 

0^  for  large  R  are  described  by 


"x  =  ike  ;  AVTS 


2  1 

°X  “  2R6  ;  DTTL 


2  5 

0,  =  T7TTT  »  DSL 


I5r! 


Hence,  in  the  linear  region,  the  absolute  value  type  bit  synchronizer 
offers  a  3dB  bit  sync  jitter  advantage  over  the  DTTL  and  a  0.  97  dB 
advantage  over  the  DSL  type  bit  synchronizer.  Furthermore,  the  AVTS 
appears  to  be  uniformly  better  than  the  other  two  configurations  at  all 
values  of  the  signal-to-noise  ratio.  If,  however,  one  is  allowed  the 
freedom  of  optimizing  the  window  width  in  the  DTTL,  then  the  bit  sync 
jitter  performance  of  the  AVTS  is  inferior  to  that  of  the  DTTL.  Fig.  VII-11 
illustrates  the  performance,  optimized  over  the  window  in  the  error  channel, 
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Figure  VII- 10.  A  Comparison  of  the  Performance  of  Three  Hit 
Synchronizer  Configurations  in  Terms  of  Their  Normalized  Hit 
Syne  Jitter.  (Reference  IS.)* 


•Reprinted  by  permission  of  the  Institute  of  Electrical  and  Electronics  Engineers.  Inc.,  345  East  47th  Street,  New  York.N.Y. 
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of  the  DTTL  versus  ite  corresponding  performance  when  the  window 
is  full,  i.  e.,  =  1.  This  figure,  when  compared  with  VII -11,  confirms 

the  above  statement. 

Finally,  a  comparison  of  performance  on  the  basis  of  mean  time 
to  first  loss  of  bit  sync  can  also  be  made  for  the  three  BSS  configuration, 
i.  e. ,  DTTL,  AVTS,  DSL.  As  one  might  suspect,  the  relative  rankings 
of  these  three  BSS  remain  the  same  as  in  the  variance  of  bit  sync  jitter 
comparison  just  made. 

VII-2.  Conclusions 

Several  conclusions  can  be  drawn  from  the  results  of  the  study 
thus  far. 

(1)  The  :,window-optimized"  DTTL  is  superior  to  the  AVTS  loop 
under  the  same  input  conditions.  This  is  true  when  one  considers  the 
variance  of  the  bit  sync  jitter  as  the  measure  of  performance. 

(2)  If  one  considers  bit  slippage  rate  5  as  a  measure  of  performance 
then  on  the  basis  of  the  phase  detector  characteristic  which  the  DTTL 
implements  versus  that  of  the  AVTS,  the  DTTL  should  have  a  lower  bit 
slippage  rate  if  the  window  in  the  error  channel  is  optimized  against  the 
noise. 

(3)  No  attempt  has  been  made  to  evaluate  the  acquisition  time  or 
acquisition  range  of  the  DTTL  or  AVTS  type  BSS.  This  is  something 
that  would  have  to  be  accomplished  by  computer  simulation.  However, 
in  no  noise,  one  can  conjecture  from  the  theory  given  in  Section  VI-1 
that  the  DTTL  will  perform  superior  to  the  AVTS  when  random  NRZ  data 
represents  the  input. 

In  summary,  then,  it  would  appear  that  on  the  basis  of  this  study 
that  a  properly  designed  DTTL  can  outperform  a  AVTS  in  both  the 
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tracking  and  acquisition  mode.  In  fact,  a  digital  version  of  the  DTTL  is 
recommended  as  the  candidate  system  to  construct  as  a  computer  software 
model.  Re  llts  of  computer  simulation  could  be  used  to  test  the  validity 
of  the  theory  as  well  as  serve  as  a  basis  for  comparison  of  the  performance 
of  other  BSS  of  interest.  The  digital  version  of  the  DTTL  type  BSS  has 
other  advantages: 

(1)  Data  Processing  in  Nonreal  Time 

(2)  Ease  of  Mechanization  in  the  Digital  Domain 

(3)  No  variation  in  analog  circuit  components 

(4)  Loop  Parameters  can  be  switched  during  the  acquisition  and 
tracking  mode 

(5)  Better  stability 

(6)  Reliability,  maintainability,  and  standardization 

(7)  Software  simulations  can  be  used  to  characterize  and  validate 
performance. 

The  major  disadvantage  of  constructing  the  DTTL  over  the  AVTS  appears 
to  be  that  it  is  more  complex  to  implement;  hence,  it  is  more  expensive 
to  build.  Finally,  it  should  be  noted  that  where  a  large  quantity  of  BSSs 
are  required  the  cost  of  implementing  the  DTTL  may  be  reasonable. 

Computer  software  programming  for  simulation  of  the  DTTL  is  best 
begun  by  considering  the  block  diagram  of  the  system.  The  details  of 
the  programming  can  proceed  on  the  basis  of  developing  the  behavior  of 
the  loop  in  terms  of  the  random  variables  which  appear  at  the  output  of 
the  in-phase  and  quadrature  integrate  and  dump  circuits.  Parameters 
v/hich  characterize  loop  performance  will  be  acquisition  time,  acquisition 
range,  bit  sync  jitter,  BEP  and  bit  slippage  rate. 
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VHI.  OUltliT  "OISE  IN  FM  DISCRIMINATORS 

Before  attempting  to  characterise  the  output  noise  of  an  FM  discrimin¬ 
ator  it  is  worthwhile  to  review  the  principles  of  operation  of  the  discrimin¬ 
ator  type  of  FM  receiver. 

Figure  VIII- 1  shows  the  block  diagram  of  the  complete  receiver.  In 
effect,  the  FM  discriminator  receiver  measures  the  rate  of  the  change  with 
respect  to  time,  of  the  angle  9  associated  with  the  signal  plus  noise  vector 
[(R,9)  in  polar  coordinates,  (c,  s)in  rectangular  coordinates^  as  shown  in 
Fig.  VIII-2.  The  equations  for  the  basic  voltages  in  the  FM  discriminator 
are  given  below  (see  Fig.  VIII- 1): 

1.  Input  to  the  receiver 

v  =  A  cos  t  +  J  m^  +  n(t)  (1) 

J  m  =  J  m(T)dT  (2) 

n(t)  =  Ax(t)  cos  uj  t  -  Ay(t)  sin  UJ  t  (3) 

c  c 

x  and  y  are  normal  (0,  <?)  Gaussian  r.  v.  's. 

2.  Input  to  Limiter 


=  Ac  cos 

uJt-Assin^t 

(4) 

c  c 

A  •  f 
s  =  sin  j 

m  +  y/ A 

(5) 

A 

c  =  cos 

f  m  +  x/A 

(6) 

x,  y  are  noise  with  bandwidth  uj  . 
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3.  Input  to  the  discriminator 


c  cos  t  -  s  sin  <*>  t 
c  c 


Ji  V 

Vs  +  c 


(7) 


4.  Input  to  the  low  pass  filter 

d  -lsv  sc  -  cs 

C3  =  dt  (Un  P  =  ~ 2 - 2  (8) 

8  +  C 

For  example,  for  sine  wave  modulation,  the  receiver  is 

v(t)  =  S(t)  +  n(t)  =  A  cos(uu  t  +  P  sin  tr  t)  +  Ax  cos  tu  t  -  Ay  sin  <*>  t 

c  m  c  c 


S(t)  +  n(t) 


R(t)  cos 


«■'  t  +  P  sin  uu  t  +  tan 
i  c:  m 


(y  cos  cp  -  x  sin  cp 
1  +  xcostp  +  y  sinCp 


! 


where  cp  =  P  sin  uu  t  and 
m 


/  ?  2 

R(t)  =  r(\  +  x  cos  Cp+  y  sin  +  (y  cos  cp  -  x  sin  cp ) 


The  output  is  given  by 


d_ 

dt 


P  sin  t  +  tan 
m 


It 


y  cos  cp  -  x  sin  cp 
+  x  cos  cp  +  y  sin  cp, 


PUJ  cos  UJ  t  +  —  Ty  cos  fP  -  x  sin  cpl  +  spikes  +  doublets 
m  m  dt  7 


n  =  gaussian  noise 


with  the  first  two  terms  being  the  signal  and  the  Gaussian  noise  respectively. 
In  what  follows  we  shall  investigate  the  statistical  characteristics  of  the 
spike  component  present  in  the  output  as  well  as  discuss  the  qualitative 
apsects  of  the  outputs. 
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Fig.  VHI-3  shows  typical  input  and  output  waveforms  for  an  FM 

discriminator  at  high  and  low  P  with  the  noise  bandwidth  (u>  )  and  the 

n 

carrier-to-noise  ratio  (CNR)  fixed.  In  addition,  the  modulation  shown  is 

assumed  to  be  sinusoidal,  i.  e. ,  m  =  Puu  cos  m  t  =  cos  U)  t  since  P  - 

m  m  n  m 

u>  /uu  .  Studies  of  these  waveforms  (which  are  readily  observed  on  an 
n  m 

oscilloscope)  have  led  various  investigators  such  as  Rice  (Ref.  19),  Cohn 
(F.ef.  20/  and  Hill  (Ref.  21)  to  describe  the  output  of  the  FM  discriminator 
as  consisting  of  the  following  components: 

1 .  Modulation 

2.  Gaussian  Noise 

3.  Doublets  (with  area  zero)  or  Additive  Noise 

4.  Spikes  (with  area  2TT)  or  Subtractive  Noise  (sometimes  called  clicks) 
The  first  component  is  the  transmitted  information.  The  last  three 

components  are  the  three  types  of  noise  which  can  occur  with  various  signal 
vector-noise  vector  configurations.  Fig.  VHI-4  shows  the  type  of  situation 
which  produces  Gaussian  noise  in  8  and  8,  i.  e. ,  when  the  amplitude  of  the 
ncise  vector  is  much  less  than  the  amplitude  of  the  signal  vector  (fixed  at 
unity  amplitude).  For  simplicity,  the  modulation  has  been  set  to  zero. 

When  the  noise  vector  is  large  and  comes  close  to  the  origin  (0),  but  does 
not  circle  the  origin  as  in  Fig.  VIII- 5,  the  discriminator  output  contains  a 
doublet.  Note  that  the  net  change  in  8,  8^  -  8  ,  is  <^217.  When  the  noise 
vector  is  large  and  does  circle  the  origin,  we  have  a  spike.  This  condition 
is  shown  in  Fig.  VIII-6.  Note  that  in  this  case  the  net  change  in  6  is 
approximately  2rr. 
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The  approximate  spectrums  of  each  of  the  four  components  at  the  output 
of  the  discriminator  are  shewn  in  Figs.  VIII-7a,  b,  c  and  d.  The  spectrum  of 
of  the  modulation  is  chosen  as  rectangular  for  simplicity.  The  spectrum  of  the 
Gaussian  noise  is  obtained  by  noting  that  the  Gaussian  noise  is  approximately 


n 

g 


(9) 


at  the  output  of  the  discriminator.  Evaluating  the  autocorrelation  function 
of  y  cos  J  m  -  x  sin  J  m,  we  find  that 


(10a) 


=  2R  (t)R  r 
x  cosjm 
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since 


E(x1x2)  =  E(yly2)  =  RX(T) 
Efx^)  =  E(x2y1)  =  0 
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m. 


For  simplicity,  one  can  assume  that 


and  that 
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The  latter  assumption  is  based  on  the  fact  that 


and  that  the  bandwidth  of  the  modulated  carrier  is  the 
of  the  noise.  Thus, 
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where  the  111  factor  is  introduced  by  the  d/dt  operator 
the  convolution  in  (17),  one  finds  that 
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(12) 

r  r  (t) 
cosj  m  (13) 
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same  as  the  bandwidth 

(17) 

in  (9).  Evaluating 
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(18) 
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The  spectra  of  the  spike  noise  and  the  doublet  noise  are  not  known 
exactly.  The  best  one  can  do  is  approximate  them  by  assuming  that  the 
spikes  and  doublets  have  typical  shapes  like  those  shown  in  Fig.  YIII-8. 
[However,  wide  variations  from  these  shapes  occur.]  The  doublet  waveform 
indicates  that  the  spectrum  is  bandlimited  and  has  a  parabolic  shape  at  the 
origin.  Thus  its  presence  can  be  largely  removed  by  the  output  filter. 

The  spike  waveform  indicates  that  the  spectrum  is  also  bandlimited  but 
is  flat  at  the  origin.  It  is  assumed  in  both  cases  that  the  sequences  of 


Figure  V III— S.  Typical  Shapes  of  Spikes  anil  Doublets. 
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doublets  and  spikes  at  the  output  of  the  di~''rimin? .or  occur  randomly 
and  have  a  finite  probability  of  occurrence.  Hence,  their  occurrence 
is  characterized  by  a  probability  density  function. 

If  we  let  N(T)  be  a  random  variable  which  denotes  the  number  of  spikes 
that  occur  in  time  T.  Assuming,  (in  accordance  with  Rice)  that  the  spikes 
are  independent  and  that  those  tending  to  increase  (decrease)  9  by  2rr  form 
a  Poisson  process  with  rate  of  occurrence  N+(N  ).  In  general,  when 
modulation  is  present,  and  N  are  not  equal.  The  probability  density 
of  z  =  N(T)  is  given  by 

OB 

p(z)  =  exp[-(N+  +  N  )T ]  6(z-k) 

k=-“ 


N  \ 

ir )  y2T'®7T'  <191 


as  can  be  shown  by  forming  the  discrete  convolution  of  the  densities  of  the 
positive  and  negative  spikes.  In  (19)  &(x)  is  the  Dirac  delta  function  and  1^0-0 
is  the  modified  Bessel  function  of  integer  order  k. 

If  one  is  concerned  with  the  type  of  modulation  in  which  the  instantaneous 
frequency  deviates  from  the  carrier  frequency  by  U)  rad/sec  for  a  time  T 
then  Rice  shows  that  the  average  rates  N+  and  N  when  the  noise  at  the  receiver 
input  is  Gaussian  (Ref.  19) 
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+  f  [  1  -  erf(vP  +  Pf 
m  m 


-  f  exp(- p )f  1  -  erf(f  /P/r)]} 
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and 


N  =  N+  +  f^  exp(-P) 


(21) 
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where 


Thus,  we  see  that  for  large  P  an  ever  majority  of  spikes  occur  in  the 
"negative  direction"  (f^  >0)  giving  rise  to  the  substractive  noise.  In  fact, 
for  large  P 


N  ~  0 

+ 

N  ~f  exp(-P) 
m 


f  >  0 
m 


(24) 


For  f  <  0  the  situation  is  reversed  of  course.  Thus  we  observe  that  the 
m 

effect  of  the  spikes  on  a  modulated  carrier  is  to  tend  to  make  the  measured 


frequencies  appear  "closer"  to  the  carrier  frequency  than  the  transmitted 
frequencies.  That  is,  confining  oneself,  for  the  moment,  to  only  error 
caused  by  spikes,  frequencies  transmitted  higher  (lower)  than  the  carrier 
will  be  measured  (by  the  discriminator)  to  be  at  that  frequency  or  a  lower 
(higher)  one  when  the  noise  is  small. 


If  one  is  willing  to  make  the  approximation  given  in  (24)  then  the 
probability  density  function  for  z  =  N(T)  can  be  approximated  by  the  Poisson 


distribution 

expf-N  T](N  T)N(T) 

p[N(T)] - TSffnf -  (fm>0) 

N(T)  =  0,1,2,...  (25) 


when  (f  >  0). 
m 

Next  we  characterize  the  spectral  properties  of  the  output  spike  noise. 
This  can  be  accomplished  as  follows:  Let  the  spikes  be  modeled  as  random 
pulses  of  wid*h  e,  height  ±2rr/e  and  rate  N  (spikes/second)  =  N  +  N  ,  as  in 
Fig.  VIII- 9.  The  autocorrelation  function  for  the  spike  process  is  illustrated 
in  Fig.  VIII- 10.  The  power  spectral  density  of  the  spike  process  is  given  by 
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l; iipire  V III  — ‘i.  \pprox  dilution  of  Spikes  With  ,i  Series  of  Unndnm  I’lilses. 
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I-‘ i mirt.-  \  Ill-Ill.  Viiloeorrclalion  l  imclioii  of  u  Series  of  Random  Pulses. 


If  z  is  small  then  the  output  filter  sees  a  power  spectrum,  due  to  the  spike 
noise,  on  the  order  of 

S(uu)  rj  16tt2N  (27) 

s 

so  that  the  output  power  due  to  the  spike  noise  is 

N  16tt2(N  +  N  )  f  (28) 

spike  +  m 

where  is  the  output  bandwidth.  What  this  says  is  that  the  energy  in  the 
spike  noise  is  spread  over  a  wide  frequency  range  and  this  range  is  wider 
than  the  modulation  bandwidth  and  wider  than  the  input  noise  bandwidth. 

In  fact,  the  equivalent  effect  of  spikes  on  the  performance  of  any  bit 
synchronization  system  is  to  drastically  raise  the  value  of  the  noise  spectral 
density  at  the  origin.  It  would  be  interesting  to  characterize  the  probability 
density  function  of  the  bit  sync  error  when  the  system  input  consists  of 
Gaussian  noise  plus  spikes. 
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Performance  Comparisons  of  FM  Discriminators  and  Matched  Filters 


If  one  assumes  perfect  bit  sync  the  BEP  using  the  FM  discriminator 
and  a  matched  filter  type  detector  are  compared  in  Fig.  VIII- 11  when  on 
is  transmitting  an  alternating  sequence  of  frequencies.  Notice  that  in  the 
region  of  small  i.  e. ,  P  s  0.  734)  the  matched  filter  and  the  FM  discrimin¬ 
ator  display  comparable  BEP.  This  is  in  the  region  where  the  additive 
Gaussian  noise  effects  the  BEP.  When  the  BEP  using  the  FM  discriminator 
result  ^  principally  from  spikes  0.734  <  P  <  4.24,  the  matched  filter  is  seen 


(•  it'iire  Vlll-I  I .  ('omp;iri.M>ii  ol  the  Muteherl  Filler  ;iml 
I'M  Diserimintilor.  (Reference  22.) 


•Reprinted  by  permission  of  the  Institute  of  Ilcctrical  and  Meet  run  ies  f.npincers.  Inc..  345  Past  47lh  Street,  New  York.  N.Y. 
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to  have  a  substantially  lower  BEP.  It  is  also  observed  that  the  minimum 
BEP  occurs  when  P  w  2  for  the  matched  filter  detector  while  the  minimum 
for  the  FM  discriminator  occurs  when  3*1.  Other  theoretical  and  experi¬ 
mental  comparisons  of  BEP  vs  modulation  index  are  illustrated  in  Fig.  VIII- 12 
when  bit  sync  is  perfect.  Lowering  of  the  FM  threshold  can  be  accomplished 
through 

1.  Employing  FM  with  Feedback 

a.  Discriminator  (FMFB) 

b.  PLL 

2.  Matched  filter  Detect  rather  than  Lowpass  Filter  and  Sample 

3.  Spike  Detection  and  Compensation. 

Spike  noise  detection  and  elimination  offers  an  advantage  in  flexibility  since 
it  is  not  a  demodulator,  but,  rather  a  device  which  can  be  implemented  at 
the  output  of  any  FM  discriminator.  The  improvement  in  threshold  perform¬ 
ance  seems  to  be  on  the  order  of  1.  5  to  2  dB  while  the  PLL  type  discrimin¬ 
ator  offers  approximately  a  3  dB  threshold  advantage  over  the  frequency 
discriminator.  The  combined  effects  of  spike  detection  and  elimination  with 
PLL  discriminator  would  then  lower  the  threshold  by  4  to  5  dB.  Spike 
detection  and  elimination  circuits  combined  with  FMFB  probably  offers  a 
3  to  4  dB  reduction  in  the  threshold  effect. 

Computer  Simulation 

Since  the  output  noise  in  an  FM  discriminator  is  impossible  to  char¬ 
acterize  in  any  exact  sense,  and  since  the  effects  due  to  the  noise  are 
unique  to  the  system  with  which  the  bit  synchronizer  electronics  is  to 
operate;  they  should  be  simulated  as  closely  as  is  practical  in  order  to 
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determine  the  performance  of  the  bit  synchronizer  as  well  as  bit  detector 
in  the  presence  of  such  phenomena.  A  simulation  in  which  the  signal  is 
summed  with  the  outputs  of  a  Gaussian  and  spike  (Poisson)  noise  generator 
which  are  calibrated  by  the  above  equations  may  perhaps  be  too  Coarse  a 
representation  to  determine  overall  performance;  therefore,  one  may  be 
forced  to  simulate  all  subsystems  in  the  overall  receiver  and  determine 
their  combined  performance. 

Time  domain  computer  simulation  techniques  are  directly  applicable  to 
problems  wherein  system  analysis  is  formidable.  The  analog  computer  can 
be  used  for  certain  problems  as  can  the  digital  computer.  However,  each 
method  has  its  particular  problems,  commonly  including  run  time  and 
complexity  of  problem  set-up.  Simulation  of  communication  links  employing 
FM  can  be  readily  accomplished.  Computation  can  be  done  at  baseband 
frequencies  so  that  the  sample  frequencies  are  not  determined  by  the  carrier 
frequency.  Nonsymmetric  bandpass  filters  can  be  accurately  represented  by 
using  bilinear  z-transform  methods  while  the  results  given  in  Ref.  ZO  can 
be  used  to  model  and  simulate  the  discriminator.  If  such  an  approach  i' 
used  one  avoids  having  to  model  the  discriminator  output  noise  directly 
so  that  the  thermal  noise  in  the  channel  is  all  that  is  required  to  model. 

There  exists  Gaussian  random  number  generators  which  readily  accomplish 
this  function.  Since  the  bit  synchronizers  studied  in  this  report  operate  in 
a  digital  fashion  it  appears  that  computer  simulation  of  the  bit  synchronizer 
chosen  will  be  required  to  determine  overall  performance  and  the  approach 
of  overall  receiver  simulation  will  be  a  fruitful  one.  Furthermore,  affects 
due  to  multipath  can  be  cranked  into  the  model  at  a  later  time.  Since  this 
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effort  may  be  costly  it  is  important  to  decide  whether  other  efforts  to 
improve  performance  are  not  worthy  of  considering  prior  to  modeling  and 
simulating  a  system  which  will  utlimately  undergo  change. 

In  developing  sample  data  simulation  techniques  for  the  system  it  is 
recommended  that  the  software  be  developed  into  three  groups  of  subroutines. 
Group  I  subroutines  would  be  those  required  to  set  up  the  simulation,  Group 
II  subroutines  actually  perform  a  function  within  the  simulation  and  each 
subroutine  will  compute  one  value  in  its  output  waveform  for  each  input 
sample  point.  Group  III  subroutines  are  those  used  by  Group  I. 

The  link  may  be  represented  by  the  following  block  diagram. 


OUTPUT 


riijurc  V  ill-11.  HI  1 .ink. 


The  bilinear  z  transform  can  be  used  to  develop  sample  data  subroutine  s 
for  making  possible  the  realization  of  filters  such  as  bandpass,  lowpass 
symmetric  and  unsymmetric.  The  dependence  of  the  sample  data  simulation 
on  the  RF  carrier  frequency  can  be  eliminated  by  using  analytic  signal  theory. 
This  approach  is  important  because  the  carrier  frequency  will  be  many 
times  the  highest  baseband  frequency  resulting  in  a  high  sample  rate  and 
long  computer  run  times  if  not  eliminated. 
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